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A. 


INTRODUCTION 


Contract  N000I A-76-C-05IA  and  its  predecessor,  contract  N00G1A- 
67-A-OI 12-0023,  have  provided  research  support  for  Pijfessor  Joseph 
Goodman  and  his  students  for  a period  of  ten  years.  The  research  performed 
under  the  contracts  has  spanned  the  fields  of  image  formation,  optical 
information  processing,  and  holography.  \ 

In  this  final  report  on  the  contracts, we  first  summarize  the 
research  accomplishments  of  the  past  contract  year  (section  II). 

Attention  is  then  turned  to  an  overview  of  the  accomplishments  of  the 
entire  ten  years  of  support  (section  III).  Included  are  short  discussions 
of  the  general  subjects  investigated,  a list  of  the  students  supported 
by  the  contractsand  their  current  professional  positions,  and  finally 
a list  of  the  publications  reporting  on  work  supported  in  part  or  in 
entirety  by  the  contracts. 
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ACCOMPLISHMENTS  OF  THE  PAST  YEAR 


ijor  emphasis  during  the  past  year  has  been  on  the  development 
of  new  methods  for  extremely  fast,  fully  parallel,  incoherent  optical 
data  processing  systems.  While  most  coherent  optical  data  processors 
require  information  to  be  input  onto  a real-time  (electronic  or  electro- 
optic) transparency  in  a serial  or  time-sequential  fashion,  incoherent 
processors  offer  the  possibility  of  fully  parallel  data  input  at 
enormous  speeds.  ^ 

The  particular  approach  to  incoherent  optical  processing  formu- 
lated and  partially  demonstrated  by  us  during  the  past  year  is  summari- 
zed in  the  Appendix,  which  is  a pre-print  of  an  article  to  appear  in 
Optics  Letters  in  January  1978.  We  believe  that  the  utilization  of 
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modern  opto-electron ic  technology  (including  LED's,  laser  diodes, 
fibers,  and  fast  photodetectors)  in  optical  signal  processors  holds 
enormous  potential  for  the  future,  and  we  hope  that  our  initial  work 
in  this  area,  as  carried  out  under  ONR  sponsorship  during  the  last  year, 
will  lead  to  a new  family  of  optical  processors  with  unique  and  useful 
properties. 

111.  OVERVIEW  OF  THE  ACCOMPLISHMENTS  OF  THE  PAST  TEN  YEARS 
A.  Areas  of  Research 

(1)  Holographic  Imaging  through  Atmospheric  Turbulence. 

The  first  area  of  research  supported  by  this  contract 

concerned  the  effects  of  atmospheric  turbulence  on  holograms  and 
holographic  images  formed  over  long  propagation  paths.  The 
research  consisted  of  both  analytical  studies,  using  the 
Tatarski  theory  of  propagation  through  turbulence,  and  experi- 
mental studies  of  the  effects  of  the  atmosphere  on  fringe 
formation  over  long  outdoor  paths.  The  research  culminated  in 
two  publications  in  the  Journal  of  the  Optical  Society  of  America 
One  graduate  student  was  supported  for  this  work,  and  he  received 
his  Ph.D.  degree  with  his  thesis  on  this  subject  matter. 

(2)  Digital  Image  Formation  from  Holographic  Data 

In  1968  we  undertook  a research  project  on  the  formation 
of  images  from  detected  holographic  data  by  means  of  digital 
computation.  The  prime  area  of  application  was  anticipated  to  be 
acoustic  holography,  for  which  a decision  between  optical  image 
formation  and  digital  image  formation  was  a very  real  one.  We 
were  successful  in  producing  the  first  digitally  reconstructed 


- 3 - 


I 

I 

E 


images  from  real  holographic  data.  In  addition,  a theoretical 
study  led  to  the  first  analytical  understanding  of  the  effects 
of  Fourier -domain  quantization  on  holographic  images.  Three 
publications  resulted,  and  one  student  received  the  Ph.D.  degree 
based  on  this  work. 

(3)  High  Efficiency  Volume  Holograms 

During  1968  and  1969.  we  undertook  a theoretical  study  of 
the  efficiency  with  which  thick  holograms  transfer  light  to  their 
reconstructed  images.  While  a theory  published  earlier  had  dealt 
with  the  theory  of  thick  holograms  consisting  of  a single 
sinusoidal  component,  our  studies  concerned  the  case  of  multiple 
stored  sinusoidal  gratings.  A paper  on  this  subject  was  presented 
at  a meeting  of  the  Optical  Society  of  America,  and  a student 
received  the  Ph.D.  degree  based  on  his  thesis  in  this  area. 

(i»)  Pre-Detection  Image  Processing  for  Enhanced  Image 

Restoration 

In  1970-1971  we  explored  a number  of  possible  ways  to  modify 
the  pupil  of  an  optical  system  in  order  to  make  the  detected 
images  more  amenable  to  post-detection  image  enhancement.  This 
work  pointed  out  for  the  first  time  that,  if  images  are  expected 
which  have  been  subjected  to  a known  degradation  (e.g.,  a focusing 
error),  a proper  apod izat ion  of  the  imaging  pupil  can  result 
in  a detected  image  which  is  especially  well  suited  to  post- 
detection  restoration.  The  work  resulted  in  a publication  and 
one  student  received  his  Ph.D.  based  on  this  work. 
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(5)  Low- Light -Level  Limitations  in  Fringe  Measurement 


and  in  Image  Restoration 

In  1970  and  1971.  a study  of  the  photon-limitations  to 
fringe-parameter  measurement  was  undertaken.  The  work  was 
motivated  by  the  interest  that  then  existed  in  a variety  of  optical 
array  techniques  for  image  enhancement,  all  of  which  required 
the  accurate  measurement  of  fringe  parameters  at  low  light  levels. 

A comprehensive  theory  was  successfully  developed.  In  addition, 
the  study  included  an  examination  of  the  photon-noise  limitations 
of  image  enhancement.  The  results  of  these  studies  were  published, 
and  one  student  received  his  Ph.D.  degree  with  a thesis  in  this 
area. 

(6)  Space-Variant  Digital  Restoration  of  Images 

In  1972,  a study  of  the  types  of  image  degradations  intro- 
duced by  relative  motions  of  a camera  and  the  subject  being 
photographed  was  undertaken.  Rather  complicated  space-variant 
blurs  were  found  to  exist.  A new  method  for  restoring  such 
blurred  images  was  developed  and  demonstrated.  Since  the  blurs 
encountered  are  space-variant,  straightforward  Fourier  domain 
restoration  procedures  could  not  be  applied,  and  more  sophisti- 
cated approaches  were  found  to  be  necessary.  The  work  resulted 
in  several  publications  and  a Ph.D.  thesis. 

(7)  Restoration  of  Degraded  Images  Formed  in  Partially 

Coherent  Light. 

The  problem  of  restoring  degraded  images  formed  in  partially 
coherent  light  is  an  extremely  difficult  one,  due  to  the  non- 
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linear  nature  of  the  image  forming  process  in  this  case.  In 
1973  we  developed  a new  approach  to  this  problem  and  demonstrated 
its  feasibility  experimentally.  A considerable  amount  of  new 
theoretical  understanding  of  the  image  forming  process  was 
developed  during  the  course  of  this  work.  The  results  were 
published,  and  one  student  received  his  Ph.D.  based  on  his 
thesis  in  this  area. 

(8)  Computer  Generated  Spatial  Filters  with  High  Light 

Efficiency 

During  1973  and  197**.  we  developed  several  new  methods  for 
constructing  computer-generated  spatial  filters  with  high  light 
efficiency.  Several  of  these  methods  were  demonstrated  experi- 
mentally. A Ph.D.  thesis  on  this  subject  resulted. 

(9)  Laser  Speckle 

During  1975  we  published  two  new  theoretical  results 
pertaining  to  laser  speckle.  One  was  the  first  analytical  deri- 
vation of  the  intensity  statistics  that  result  from  the  addition 
of  N partially  correlated  speckle  patterns.  The  second  was  an 
analytical  derivation  of  the  intensity  statistics  obtained  in  a 
speckle  pattern  generated  by  a surface  which  has  a roughness 
that  is  smaller  than  a wavelength.  Following  these  publications, 
we  prepared  a review  article  on  the  statistics  of  laser  speckle 
patterns.  This  article  is  now  one  of  the  standard  references  in 
the  field. 

(10)  Nonlinear  Optical  Data  Processing 

During  the  years  1973~1976  we  developed  new  methods  for 
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performing  nonlinear  optical  data  processing  operations  using 
coherent  optical  systems.  These  methods  are  based  on  the  use 
of  half-tone-screen  encoding  processes.  The  work  had  sufficient 
impact  to  motivate  a number  of  other  research  groups  in  the 
U.S.  and  abroad  to  continue  efforts  along  these  lines  based  on 
our  early  publications.  Our  work  included  methods  for  realiz- 
ing logarithmic  transformations,  analog-to-d ig i tal  conversion, 
and  pseudo-color  picture  encoding.  Approximately  six  publi- 
cations resulted  from  this  work. 

(11)  Incoherent  Matrix-Vector  Multiplier 

This  work  has  begun  in  1977,  and  is  reported  on  in  detail 
in  the  Append ix. 

In  addition  to  the  above  research  areas  specifically  funded  by  our 
ONR  contracts,  several  other  research  areas  were  initiated  with  the  help 
of  this  contract,  and  then  funded  for  more  in-depth  research  by  other 
agencies.  These  areas  include: 

Interferometric  imaging  methods  for  astronomical  imaging 
through  turbulence  (NSF); 

Computer  generated  holographic  memories  (NSF);  and 
Space-variant  optical  data  processing  systems  (NSF). 

Four  students  received  their  Ph.O's  under  these  latter  research  efforts. 

B.  Students  Receiving  Ph.O's  with  Partial  or  Total  Support  from 
our  ONR  Contracts 

I.  NAME:  Jack  0.  Gaskill  (total  support) 

Thesis  title  and  date:  "Holographic  Imaging  Through  a 


Randomly  Inhomogeneous  Medium",  May  1 968 - 

Present  position:  Administrator  for  Academic  Affairs  and 

Professor  of  Optical  Sciences,  Optical  Sciences  Center, 

University  of  Arizona,  Tucson,  Arizona. 

NAME:  Antonio  Silvestri  (total  support) 

Thesis  title  and  date"  "Digital  Image  Formation  and  Fourier 
Domain  Quantization",  May  1971. 

Present  position:  Project  Leader,  Electromagnetic  Systems 
Laboratories,  Inc.,  Sunnyvale,  Ca. 

NAME:  Richard  A.  Baugh  (partial  support) 

Thesis  title  and  date:  "High-Efficiency  Volume  Holography", 

June  1969- 

Present  position:  Staff  Scientist,  Hewlett-Packard  Co., 
Cupertino,  Ca. 

NAME:  F.  Donald  Russell  (partial  support) 

Thesis  title  and  date:  "Predetection  and  Postdetection  Filtering 
for  Improved  Resolution  in  Optical  Systems",  August  1970. 

Present  position:  Manager,  Ford  Aeronutronic,  Palo  Alto,  Ca. 
NAME:  John  F.  Walkup  (total  support) 

Thesis  title  and  date:  "Limitations  in  Interferometric 
Measurements  and  Image  Restoration  at  Low  Light  Levels", 

July  1971- 

Present  position:  Associate  Professor  of  Electrical  Engineering, 
Texas  Tech  University,  Lubbock,  Texas. 
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NAME:  A. A.  Sawchuk  (partial  support) 

Thesis  title  and  date:  "Space-Variant  Image  Motion  Degradation", 
June  1972. 

Present  position:  Associate  Professor  of  Electrical  Engineering, 
University  of  Southern  California,  Los  Angeles,  Ca. 

7.  NAME:  Mete  Severcan  (partial  support) 

Thesis  title  and  date:  "Computer  Generation  of  Coherent  Optical 
Filters  with  High  Light  Efficiency  and  Large  Dynamic  Range", 
Present  position:  Instructor,  Middle  East  Technical  University, 
Ankara,  Turkey. 

8.  NAME:  Kalyan  Dutta  (total  support) 

Thesis  title  and  date:  "Sampling  and  Restoration  of  Images 
Formed  in  Partially  Coherent  Light",  December  1 97^* - 
Present  position:  Staff  Scientist,  Block  Engineering, 

Cambridge,  Mass. 


C . Publications  on  Work  Supported  in  Part  or  in  Whole  by  our 
ONR  Contracts 


1.  J.W. Goodman,  "Digital  Image  Formation  from  Electronically 
Detected  Holograms",  Proc.  Seminar  on  Computerized  Imaging 
Techniques,  Soc.  of  Photo-Optical  Instrumentation  Engineers, 
pp . XV 1-1  through  XVI-6,  1967  (Monograph). 

2.  J.W. Goodman  and  R.W.  Lawrence,  "Digital  Image  Formation  from 
Electronically  Detected  Holograms",  Appl .Phys.Lett. , 1 1 , 

pp. 77-79,  August  1967. 

3.  J.W. Goodman,  "Systems  Applications  of  Holography",  Proc . 
Seminar  on  Holography,  Soc.  of  Photo-Optical  Instrumentation 
Engineers,  pp.l*»7-152,  (Monograph),  1968. 

k.  J.W. Goodman,  R.B.  Miles  and  R.B.  Kimball,  "Comparative 
Noise  Performance  of  Photographic  Emulsions  in  Conventional  and 
Holographic  Imagery",  J . Opt . Soc .Am. , 58,  pp.609-6l**,  May  I968. 
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5.  J.O.  Gaskill,  "Imaging  through  a Randomly  Inhomogeneous 
Medium  by  Wavefront  Reconstruction",  J . Opt . Soc . Am. , 58, 

pp. 600-608,  May  1 968 . 

6.  J.W. Goodman  and  A.  Silvestri,  "Digital  Reconstruction  of 
Holographic  Images",  NERtM  RECORD,  1 0,  p . 1 1 8 , 1 968 . 

7.  J.W. Goodman  and  G.R.  Knight,  "Effects  of  Film  Nonlinearities 
in  Wave-f ront-Reconstruct ion  Imaging  of  Diffuse  Objects", 

J. Opt. Soc. Am. , 58,  pp. 1 276-1 283 , September  1968. 

8.  J.W. Goodman,  "Digital  Image  Formation  from  Detected 
Holographic  Data",  Acoustical  Holography,  (ed.  A.F.  Metherel  l , 

H.M.A.  El-Sum  and  L.  Larmore),  Plenum  Press,  pp . 1 73 - 1 85 > 1 969 . 

9.  J.W. Goodman,  "Use  of  a Large  Aperture  Optical  System  as  a 
Triple  Interferometer  for  Removal  of  Atmospheric  Image  Degra- 
dations", Evaluation  of  Motion  Degraded  Images,  NASA  publication 
SP-193,  1959: 

: i 

10.  J.D.  Gaskill,  "Atmospheric  Degradation  of  Holographic 
Images",  J .Opt. Soc. Am. , 59,  pp.308~3l8,  March  1 969 * 

11.  J.D.  Gaskill  and  J.W. Goodman,  "Use  of  Multiple  Reference 
Sources  to  Increase  the  Effective  Field  of  View  in  Lensless 
Four ier-Transform  Holography",  Proc.  IEEE,  57.  p-823,  May  1969- 

12.  J.W. Goodman,  "Synthetic-Aperture  Optics",  in  Progress  in 
Opt i cs , Vol . VIII,  (ed.  E.  Wolf),  North  Holland  Publishing  Co., 
Amsterdam,  pp.1-50,  1970. 

13.  J.W. Goodman,  "Analogy  Between  Holography  and  Interferometric 
Image  Formation",  J .Opt ■ Soc .Am. , 60,  pp. 506-509,  April  1970. 

14.  J.W. Goodman  and  A.M.  Silvestri,  "Some  Effects  of  Fourier- 
Domain  Phase  Quantization",  IBM  J. Res. Dev. , 1 4 , pp. 478-484, 

September  1970. 

15.  P.D,  Russell  and  J.W. Goodman,  "Nonredundant  Arrays  and  Post- 
detection Processing  for  Aberration  Compensation  in  Incoherent 
Imaging",  J. Opt .Soc. Am. , 61 , pp . 1 82- 1 9 1 » 

16.  J.W. Goodman,  "An  Introduction  to  the  Principles  and  Appli- 
cations of  Holography",  Proc.  IEEE,  59,  pp . 1 292- 1 304 , September  1971. 

17-  A. A.  Sawchuk,  "Space-Variant  Image  Motion  Degradation  and 
Restoration",  Proc.  IEEE,  60,  pp. 854-861 , July  1972. 

18.  J.F.  Walkup  and  and  J.W.  Goodman,  "Limitations  in  Fringe 
Parameter  Estimation  at  Low  Light  Levels",  J. Opt. Soc. Am. , 63, 

PP. 399-407,  April  1973. 
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19-  J.W. Goodman,  "Statistical  Properties  of  Mutual  Intensity 
with  Finite  Measurement  Time",  Appl . Phys . 2,  pp.95-l01, 

August  19/3* 

20.  H.  Kato  and  J.W.  Goodman,  "Nonlinear  Transformations  and 
Logarithmic  Filtering  in  Coherent  Optical  Systems",  Opt ics 
Comm. , 8^  pp. 378-381,  August  1 973  - 

21.  J.W. Goodman,  "Probability  Density  Function  of  the  Sum  of 
N Partially  Correlated  Speckle  Patterns",  Optics  Comm.  13. 
pp. 244-247  (1975). 

22.  J.W. Goodman,  "Dependence  of  Image  Speckle  Contrast  on 
Surface  Roughness",  Optics  Comm.  14,  pp. 324-327  ( 1 97 5 ) • 

23.  H.  Kato  and  J.W.  Goodman,  "Nonlinear  Filtering  in  Coherent 
Optical  Systems  through  Halftone  Screen  Processes",  Appl. Opt. 

T4,  pp. 181 3-1824  (1975). 

24.  J.W. Goodman,  "Statistical  Properties  of  Laser  Speckle 
Patterns",  in  Laser  Speckle  and  Related  Phenomena,  Springer- 
Ver lag,  Heideiberg,  pp.9"75  (1975). 

25.  J.W. Goodman,  "Noise  in  Coherent  Optical  Information 
Processing",  in  Optical  Information  Processing,  Plenum  Press, 
(1976). 

26.  J.W. Goodman,  "Operations  Achievable  in  Coherent  Optical 
Data  Processing",  in  Proceedings  of  the  1975  Electro-Optical 
System  Design  Conference,  Industrial  and  Scientific  Conference 
Management,  Inc.,  Chicago,  111.  ( 1 975) - 

27.  H-K.  Liu,  J.W.  Goodman  and  J.L-H.  Chan,  "Equ idens i tometry 
by  Coherent  Optical  Filtering",  Appl .Optics,  15,  pp. 2394-2399. 
October  1976. 

28.  H-K.  Liu  and  J.W. Goodman,  "A  New  Coherent  Optical  Pseudo- 
Color  Encoder",  Nouvelle  Revue  d'Optique,  pp.285“289. 

Sept. -Oct.  1976. 

29.  J.W. Goodman,  "Operations  Achievable  with  Coherent  Optical 
Information  Processing  Systems",  (invited  paper),  Proc.  IEEE,  65., 
pp. 29-38,  January  1977. 

30.  J.W. Goodman,  "Some  Fundamental  Properties  of  Speckle", 

J.  Opt.Soc.Am.,  66,  pp. 1 145- 1 l 50,  November  1976. 

31.  H-K.  Liu  and  J.W.  Goodman,  "Pseudo-Color  Encoding  of 
Pictures  by  Means  of  Halftone  Screens  and  Coherent  Optical 
Filtering",  Proc.  SPIE,  83 , pp. 125-129,  San  Diego,  Ca.  August 
24-26,  1976. 
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32.  K.  Dutta  and  J.W.  Goodman,  "Reconstruction  of  Images  and 
Partially  Coherent  Objects  from  Samples  of  Mutual  Intensity", 
J.Opt.Soc.Am. , 67,  pp. 796"803,  June  1977. 

33.  J.W. Goodman,  "Some  Fundamental  Properties  of  Speckle", 
invited  lecture,  International  School  of  Quantum  Electronics: 
Coherent  Optical  Engineering,  Tuscany,  Italy,  September  1976, 
to  be  published  by  North  Holland  Publishing  Co. 

3**.  J.W. Goodman,  "Coherent  Optical  Image  Deblurring",  invited 
lecture,  International  School  of  Quantum  Electronics:  Coherent 
Optical  Engineering,  Tuscany,  Italy,  September  1976,  to  be 
published  by  North-Hoi  1 and  Publishing  Co. 

35-  J.W. Goodman,  A.R.  Dias,  and  L.M.  Woody,  "A  Ful ly-Paral lei , 
High-Speed  Incoherent  Optical  Method  for  Performing  Discrete 
Fourier  Transforms",  accepted  for  publication  in  Optics  Letters 


I V . CONCLUDING  REMARKS 

The  Principal  Investigator  for  this  contract,  Professor  Joseph 
W.  Goodman,  would  like  to  express  his  personal  thanks  to  the  Physics 
Branch  of  the  Office  of  Naval  Research  for  their  support  over  the  past 
ten  years.  This  contract  has  contributed  enormously  to  his  own 
professional  growth,  and  has  further  provided  partial  or  full  support 
for  eight  graduate  students  at  the  beginning  of  their  respective 
professional  careers.  We  believe  that  the  research  accomplished  under 
the  contract  has  had  major  impact  on  several  different  areas  of  physics 
and  engineering. 

The  freedom  afforded  to  the  Principal  Investigator  to  select  and 
pursue  those  research  efforts  he  perceived  to  be  of  major  importance 
is  especially  appreciated. 
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APPENDIX 


A FULLY-PARALLEL,  HIGH-SPEED 
INCOHERENT  OPTICAL  METHOD  FOR 
PERFORMING  DISCRETE  FOURIER  TRANSFORMS* 


J.W.  Goodman 
A.R.  Dias 
L.M.  Woody 


Department  of  Electrical  Engineering 
Stanford  University 
Stanford,  California  9^305 


An  incoherent  optical  data  processing  method  is  described  which 
has  the  potential  for  performing  discrete  Fourier  transforms  of  short 
length  at  rates  far  exceeding  those  afforded  by  both  special  purpose 
digital  hardware  and  representative  coherent  optical  processors. 


* 


Work  supported  by  the  Office  of  Naval  Research. 


We  report  here  on  an  incoherent  optical  method  for  performing 
discrete  Fourier  transforms  (DFT's)  which  has  the  potential  for  extremely 


high  data  throughput  rate.  The  DFT  operation  may  be  viewed  as  a process 
of  multiplying  an  input  vector  f (consisting  of  N possibly  complex- 
valued input  samples)  times  an  N*N  matrix  JJ,  (the  n.m1^1  element 
being  exp(- j2nnm/N) ) to  yield  an  output  vector  g (consisting  of  the 
N complex  Fourier  coefficients);  thus  we  desire  to  perform 

g = Ji  ? . (I) 

Two  separate  issues  must  be  addressed  in  describing  the  method  of 
interest  here:  (1)  How  do  we  perform  the  matrix  product  in  a highly 
parallel  and  fast  way?  (2)  How  do  we  perform  complex  arithmetic  using 
incoherent  light,  for  which  only  non-negative  and  real  quantities 
(intensities)  can  be  manipulated. 

To  address  the  first  issue,  suppose  that  the  elements  of  t and  Jj, 
are  non-negative  and  real.  Then  the  system  shown  in  Fig.  1 can  be  used 
to  perform  the  matrix-vector  product.  The  elements  of  T are  entered 
in  parallel  by  controlling  the  intensities  of  N light  emitting  diodes 
(LED's).  Lenses  Lj  and  image  the  LED  array  horizontally  onto  the 

matrix  mask  M , while  spreading  the  light  from  any  single  LED  vertically 
to  fill  an  entire  column  of  the  matrix  mask.  Lens  L^  is  a field  lens. 
The  matrix  mask  M consists  of  NxN  subcells,  each  containing  a trans- 
parent area  proportional  to  one  of  the  matrix  elements.  Lens  L^  is  a 
cylindrical  lenslet  array,  which  is  not  essential  to  the  operation  of  the 
system,  but  which  can  be  used  to  improve  light  efficiency.  Lens  combin- 


ation L,.  collects  all  light  from  a given  row  and  brings  it  to  focus  on 
one  element  of  a vertical  array  of  N photodetectors.  Each  photodetector 
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measures  the  value  of  one  component  of  the  output  vector  g. 

To  permit  the  multiplication  of  a matrix  with  complex  elements 

times  a vector  t with  complex  elements,  we  decompose  each  of  these 

1 2 

quantities  as  follows  ’ : 


? = ?(0)  +?(,)  exp(j 2ir/3)  + ?(2)  exp(jW3) 

H - lH,(0)  + lH((,)  exp(J2ir/3)  +JL(2)  exp(jW3) 

where  r(°)  , f(D 

and  *<2) 

consist  of  N real  and  non-negative 
elements  each,  and  ^ ancl  consist  of  NxN  real 

and  non-negative  elements.  If  the  output  vector  g is  similarly 
decomposed,  then  we  find  that  the  overall  matrix-vector  product  can 


be  expressed  as 


ji(0)  d2)  f ?<°> 

„(1)  „(0)  (2)  ,(1) 

» **  i >n_j  I 


*,<'>  T<2> 


Thus  complex  operations  can  be  performed  at  a price  of  a factor  of  three 
in  the  length  of  the  input  and  output  vectors. 

Simple  electronic  circuits  for  producing  the  components  ?<°)  , ?(,) 
and  ?(2)  from  ? exist1,  as  do  likewise  simple  circuits  for  producing 
the  real  and  imaginary  parts  of  g from  g ^ , g^  and  g^.. 

Experiments  have  been  carried  out  to  verify  the  ability  to  perform 
complex  arithmetic.  The  source  was  an  unfiltered,  linear  filament,  clear 
envelope,  incandescent  bulb.  The  30x30  matrix  mask  used  to  perform  a 
10-point  OFT  is  shown  in  Fig.  2.  This  mask  is  designed  so  that  the 
three  entire  vectors  and  are  entere£|  side  by  side, 
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whereas  the  three  output  components  and  g^  for  the 

k*^1  Fourier  coefficient  appear  side  by  side.  Thus  the  output  display 
shows  each  DFT  component  as  a triplet  of  real  and  nonnegative  components. 

For  this  experiment  the  input  functions  were  entered  by  hand  as 
masks  placed  against  the  matrix  mask,  and  output  functions  were  detected 
on  a 1021*  element  Reticon  CCD  detector  array.  Figure  3 shows  both 
theoretical  output  distributions  and  experimentally  obtained  output 
distributions,  the  latter  being  photographed  from  an  oscilloscope  display. 

In  parts  (a)  and  (b),  the  function  to  be  transformed  consists  of  the 
sequence  (1 ,0,0, 0,0, 0,0, 0,0,0) . The  resulting  DFT  should  be  entirely  real 
and  of  constant  magnitude.  As  shown  in  these  figures,  the  DFT  components 
along  the  real  axis  are  all  non-zero  and  equal,  while  the  componenents 
along  120°  and  240°  are  all  zero. 

In  parts  (c)  and  (d),  the  input  sequence  was  entirely  real  and 
constant.  The  DFT  consists  of  a large  real  zero-frequency  component  (on 
the  far  right),  followed  by  triplets  of  equal  strength  for  all  other 
DFT  components.  Some  thought  shows  that  any  DFT  component  with  elements 
g^  , g^  and  g ^ exactly  equal  is  equivalent  to  a zero  result. 

Hence  all  DFT  components,  except  the  zero  frequency  component,  are  zero. 

Parts  (e)  and  (f)  show  the  results  when  the  entire  matrix  mask  is 
uniformly  illuminated.  In  this  case  some  thought  shows  that  the  input  is 
effectively  a sequence  containing  all  zeros.  The  output  DFT  shows 
triplets  of  equal  strength,  or  a sequence  of  all  zeros  for  the  output. 

A system  composed  of  96  high  speed  LED's  and  96  ava lance  photo- 
diodes would  be  capable  of  performing  a 32-point  DFT.  Commercially 
available  components  have  sufficient  bandwidth,  output  power  and  sensitivity 
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to  permit  such  a DFT  to  be  performed  every  10  nanoseconds.  The  total 

q 

throughput  rate  for  such  a processor  is  about  3 x 1 0 complex  samples  per 
second,  while  a corresponding  number  for  special-purpose  digital  array 
processors  is  about  3*10^  complex  samples  per  second,  and  a represent- 
ative coherent  optical  processor3  has  a throughput  of  3*10^  real  samples 
per  second. 

The  chief  significance  of  this  processor  is  that  the  input  data 
can  be  entered  in  parallel,  and  it  is  this  fact  that  leads  to  its  high 
throughput  rate.  Another  system  recently  described4’5  performs  a 
similar  matrix-vector  product,  but  the  data  must  be  entered  serially, 
and  as  a consequence  the  throughput  rate  is  much  lower.  The  processor 
described  here  is  especially  well  suited  for  problems  in  which  the 
elements  of  the  input  vector  are  gathered  by  parallel  sensors. 

Of  course,  matrices  other  than  the  DFT  matrix  can  also  be  used  if  desired. 
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